Finding the remnants of lost jets at RHIC 
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A fast parton propagating through partonic matter loses energy by medium-induced gluon radi- 
ation. We propose a method to detect the energy loss of jets and the final-state interaction of the 
radiated gluons and the jets with the partonic medium. By embedding jets and radiated energy 
into a parton cascade, we find that the fate of lost jets can be studied by analyzing the momentum 
distribution of soft hadrons relative to the jet axis. 

PACS numbers: 12.38.Mh, 24.85.+p, 25.75.-q 



The attenuation of hard jets produced in pQCD pro- 
cesses provides an excellent opportunity to study the 
properties of hot and dense partonic matter expected 
to be formed in ultrarelativistic heavy ion collisions 
Q, 0i S Due to the energy lost by the jets as a re- 
sult of multiple elastic and inelastic interactions in the 
partonic matter, the final hadron px distribution from 
the fragmenting jets will be suppressed compared to that 
from nucleon-nucleon reference spectra convoluted with 
the number of binary collisions. Measurement of high px 
hadrons is thus expected to provide a diagnostic tool for 
the total energy loss AE of jets in the dense partonic 
medium 0, @ . Relative to the binary collision scaled pp 
reference spectrum, experimental data 0, 0, 0, 0] for 
central Au+Au collisions at RHIC indeed showed a sup- 
pression for 7r°s and ft. ± , where the suppression is small- 
est at px ~ 2 GeV/c and increases to an approximately 
constant value for pt — 4 — 10 GeV/c. Apart from the 
suppression of high px hadrons at RHIC, the suppres- 
sion of the away-side jets compared to the near-side jets 
at high pt 0] is an indication of jet quenching for the 
away-side jet. 

Even the disappearance of away-side jets might be con- 
sidered as only circumstantial evidence of jet-energy loss. 
For instance, one could imagine a scenario where mono- 
jets are produced due to the 2 — ► 1 gluon fusion pro- 
cesses, e.g., in the gluon saturation model 0. If the 
away-side jet were truly produced then absorbed by the 
medium, the extra momentum and energy should sur- 
vive. Our aim is to investigate the degree to which the 
lost energy and momentum can be extracted from low-pr 
observables. We propose a straight-forward strategy of 
gating on the observation of a high-energy parton, then 
analyzing the accompanying alteration of the low-energy 
spectra relative to the azimuthal direction of the jet. The 
distribution of these extra particles, as a function of px 
and as function of their rapidity and azimuthal angle rel- 
ative to the triggering jet, might provide critical insight 
for understanding the mechanism for jet-energy loss. 

To investigate these issues, we consider a baseline dy- 
namical model of the collision, onto which hard jets have 
been embedded. The base distribution consists of the HI- 



JING model 0,0] , which is used to generate a set of ini- 
tial partons, and ZPC, a parton cascade which simulates 
the subsequent elastic partonic interactions 0]. HIJING 
provides the initial momentum and space-time coordi- 
nates of partons from minijets and strings which originate 
from soft nucleon-nucleon collisions. The parton-parton 
elastic scattering in ZPC is determined by the cross sec- 
tion, 
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where t and s are the usual Mandelstam variables, 
C a = 9/8,1/2,2/9 for gg,gq,qq scatterings, and a s is 
the strong coupling constant assumed to be 0.30. The 
effective screening mass fi is assumed to independent 
of temperature and density and is used as a parame- 
ter to obtain the desired parton-parton elastic cross sec- 
tion a » Carnal/ . We have used a typical value for 
/i = 0.50 GeV. After the partons freeze out, they are 
traced back to their parent strings and are converted to 
hadrons using the Lund string fragmentation model |l6l | . 

Since the jet production cross section is small, we in- 
clude additional hard jets in order to enhance statis- 
tics. The multiplicity of the additional jets in an event 
is obtained from a Poisson distribution with a mean 
(n) = 10. The underlying dynamics of the base distri- 
bution (dJVparton/dy ~ 250 at y/s = 200A GeV) is not 
significantly affected by the smaller number of the em- 
bedded jets. This technique of oversampling has been 
widely used to study rare particle production, such as 
K~ at SIS/GSI [II, multistrange particles at SPS 0, 
and J I ifj at RHIC ll9| . In p + p collision, the lowest or- 
der (LO) pQCD momentum distribution of hard partons 
that propagate in the parton cascade is given by 
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In Eq. @ , x a , %b are the momenta fractions and kx a , 
kj 1 , are the intrinsic transverse momenta of the initial 
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scattered partons. A constant factor K = 2 is used 
to account for the NLO corrections. For the collinear 
parton distribution function (PDF) f a (a; a ,Q 2 ), we em- 
ploy the LO Gluck-Reya-Vogt (GRV94) parametrization 
[20|, and for the transverse PDF g{kx a ), we use Gaussian 
intrinsic kx distribution of width (k|, ) = 1 GeV 2 /c 2 . 
We neglect nuclear shadowing for the jets as its effect is 
rather small at mid-rapidity at RHIC energies 0, H3| • 
Since we explore events with at least one triggered hadron 
of p^ ie > 4 GeV/c, a minimum momentum transfer of 
p™ ln = 4 GeV/c for the additional jets is assumed. In the 
hard processes, a proper formation time of To — l/mx 
is used for parton production with transverse mass tut- 
The initial transverse coordinate of a hard jet is speci- 
fied by the number of binary collision distribution for two 
Woods-Saxon geometries. 

The hard jets can suffer energy loss by elastic scatter- 
ing with the partons from HIJING. However, the domi- 
nant energy loss is via medium-induced gluon radiation 
along the trajectory of the jet. We consider here the dom- 
inant first-order total energy loss expression in opacity, 
X = L/Xg, from the reaction operator approach [it E^|: 



AE = 2na d s C a C A / rfrp(r(r), t)(t - r ) log 
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Here, r(r) and E are the position and initial energy of 
the jet, p(r(r),r) is the parton density. The color fac- 
tor C a is given in Eq. Q and may be expressed as 
C a — CrCt i 'd-A, where Cr and Ct are the color Casimirs 
of the jet and the target parton, and cLa is the dimension 
in the adjoint representation with Casimir Ca- The ef- 
fective path length required by the radiated gluons for 
traversing the matter, L, is taken to be 4 fm. Because of 
jet quenching the fragmentation function (FF) for a jet c 
to a hadron h in vacuum, Dh/ C (z c , Q 2 ), is altered due to 
the modification of the kinetic variables of the jet, and 
the effective fragmentation function follows the relation 



z cD' h/c (z c , Ql) = z' c D h/c (z' c , Ql) + N g z g D h/g (zg,Q 2 g ) , 



Ph 

Pc 



Ph 



Ph 



p c - AE(p c 



AE{ Pc )/N g 



(4) 



The second term on the right eventually gives the soft 
hadrons from the fragmentation of N g radiated gluons. 
Note that the modified fragmentation function satisfies 
the sum rule J dz c z c D' h ^ c (z c ,Ql) = I- The initial hard 
jets after the partonic phase are fragmented using the LO 
Binnewies-Kniehl-Kramer (BKK) parametrization [24J of 
the fragmentation functions. 

In absence of any definitive relation for how radi- 
ated energy should be translated into particles, we take 
recourse to a statistical prescription. Assuming local 
thermodynamical equilibrium of the parton medium, if 
AE(r. At) is the energy radiated in a time interval At 
(evaluated from Eq. ©), the number of thermal gluons 



emitted in At can be related to the entropy increase AS" 
of the plasma, 
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Ng(r,Ar)^-AS=- 
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where the local temperature T w e(r, r)/3p(r, t) is ex- 
tracted from the cascade by viewing the local parton 
number and energy densities p(r, t) and e(r, t). The fac- 
tor of 1/4 comes from considering an ideal massless gas 
of classical particles in equilibrium. The radiated gluons 
are formed at space-time points along the parent jet's tra- 
jectory consistent with the energy loss described in Eq. 
© ■ They are initially directed along the axis of the par- 
ent jet and subsequently scatter via elastic collision with 
the neighboring partons after a formation time To as used 
in HIJING. Effectively, this procedure converts 100% of 
the radiated energy into thermal energy, while maintain- 
ing momentum conservation. A less-efficient prescription 
for thermalizing the radiated energy would result in a re- 
duced multiplicity. Based on local parton-hadron duality 
[25I Eif , each of these radiated gluons after freeze-out are 
converted to a pion with equal probability for the three 
charge states. 
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FIG. 1: The production rate of gluons at all rapidities from 
radiative jet energy loss in central (6 = fm) Au+Au colli- 
sions at RHIC cm. energy of ^Js = 200^4 GeV. The results 
are for jets with different pr (dashed lines) and for the total 
jets (solid line) in an event, normalized to the average num- 
ber of jets (Nj) with the given pr- The dash-dotted line gives 
the rate of radiative energy loss (dAE /dt)/ (Nj) in GeV/fm. 
The rise at small times is the result of performing the calcu- 
lation in Cartesian time, rather than proper time. The inset 
shows the energy dependence of fractional energy loss AE/E 
of gluon jets due to radiation (dashed line) and after elastic 
scattering (solid). 
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The inset of Fig. ^ shows the fractional energy loss 
AE/E for the fast gluons traversing the partonic medium 
for central Au+Au collisions at a center of mass energy 
xfs = 200A GcV. In contrast to the asymptotic BDMS 
121 energy loss for a thick static plasma, AS B dms = 
■koP s C r C a /2 -v-L?p (the factor v oc log(L/A 9 ) ~ 1 - 3), 
where AE/E increases rapidly with decreasing energy, 
the GLV 0, 113 first-order radiative energy loss AE 
(dashed line) exhibits almost a linear energy dependence 
for E = 4—7 GcV while at higher energy AE/E reveals a 
\og(E)/E decrease. Elastic collision with other partons, 
which is modeled with the parton cascade, leads to an 
additional energy loss by 25% (solid line). 

The dynamical evolution of the average energy loss per 
jet in an event (for jets at all rapidities), (dAE / dt) / (Nj) , 
rises for small times due to the longitudinal expansion of 
the system as shown in Fig. ^ If proper time had been 
used rather than Cartesian time, the rate would have 
monotonically fallen due to the falling density. 

Figure ^ also shows the rapidity integrated dynami- 
cal rate of gluon radiation per average number of jets, 
(dN g /dt)/(Nj), as evaluated from Eq. (0). The results 
are for total jets in an event, and also for jets with definite 
Pt- The emission pattern is found to follow that for the 
energy loss. Though the number of radiated gluons per 
(Nj) is seen to increase with Pt, the rapidly decreasing 
jet production cross section with increasing jet energy E 
actually leads to smallest number of gluon contribution 
from jets with px > 12 GeV/c. Although higher-energy 
jets are rare, they radiate more particles due to fact that 
they have more available energy and due to the logarith- 
mic dependence of the radiation in Eq. We find 
that partons with initial px ~ 9 GeV/c produce the bulk 
of the radiated gluons. Using the statistical prescription 
described above, about 5.5 soft gluons are emitted in an 
event per jet. 

Since jets are typically formed in pairs, the observa- 
tion of a jet can be used as a trigger. One can then 
search for the remains of the balancing jet, as well as any 
particle whose existence is derived from medium-induced 
radiation from the triggered jet. For our calculations, the 
trigger is defined as a hadron at midrapidity \y | < 1 with 
its transverse momentum p^ ls above 4 GeV/c. The x di- 
rection is then defined to be perpendicular to the beam 
axis and pointing along the direction of the triggering 
hadron. Particles radiated from the trigger jet would 
typically have positive values of p x while particles asso- 
ciated with the jet's partner would have negative values 
of p x . For a perfect detector, momentum conservation 
implies 

J Pxp(Px) dp x = -P T , (6) 

where Pt = p T lg ^ s ^ ne momentum of the trigger particle 
and p(px) = dN/dp x refers to the momentum density of 



all other particles. Of course, this constraint is signifi- 
cantly relaxed by finite acceptance. By making a similar 
distribution po from other events, without the require- 
ment of a trigger, one can define 

&p(Px) = p(Px) - Po(Px) , (7) 

which describes the variation of the momentum density 
due to the jet. Two quantities related to Ap are of par- 
ticular interest: 

AN = J Ap( Px ) d Px , 

AP X = J p x Ap{p x ) dp x . (8) 

AJV counts the net number of additional particles, while 
AP X describes the momentum, which should be a fraction 
of —Pt, determined by the acceptance. 
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FIG. 2: The distribution p(p x ) for particles produced at mid- 
rapidity \y\ < 1 in events with a triggered hadron 4 < p^ lg < 6 
GeV/c in the same rapidity range. The results are for total 
number of particles (solid circles) and for particles from ra- 
diative jet energy loss (open circles) in central (6 = fm) 
Au+Au collisions at RHIC energy of y/s = 200A GeV. The 
p + p results at the same energy are shown by open squares. 
The ratio Raa(Px) for the multiplicities in Au+Au collisions 
to that from p+p collisions for the triggered jet event is shown 
in the inset. 

Figure |21 shows p(p x ) for events with a trigger momen- 
tum, 4 < Pt = p T ie < 6 GeV/c. For large negative 
values of p x , p(p x ) is suppressed in Au+Au collisions 
relative to pp results as expected due to the suppres- 
sion of away-side jets 01 • The situation is reversed for 
\p x \ < 600 MeV/c, where the spectrum is dominated by 
the radiated and rescattered particles. The momentum 
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constraint, Eq. ©, is satisfied by an increase of low 
momentum particles to compensate for the loss of the 
away-side jet. The net number of extra particles, AN, 
is higher in the Au+Au case as would be expected from 
the induced-radiation and from the statistical nature of 
Eq. JSJ which converts the radiated energy into particles 
with momenta characterized by the temperature. The 
alteration of the jets is especially apparent in the inset 
of Fig. [3 which shows the ratio of the Au+Au and pp 
results. At large negative p x the suppression is approxi- 
mately a factor of 3, while the enhancement at low p x is 
approximately a factor of two. 

Also shown in Fig. [21 (open circles) is the distribution 
of the particles in Au+Au collisions that stems from the 
induced radiation described in Eq. J2J. As expected, 
these particles are at much lower px than those from pp 
collisions due to their thermalizing with the medium. 

For this range of triggers the number of extra parti- 
cles, both charged and neutral, produced at mid-rapidity 
was approximately AN = 19, versus 14 for pp collisions. 
We emphasize that this is a relative increase and that 
determining po is difficult, as discussed later. The num- 
bers increase by ~ 20% if all rapidities are evaluated. Of 
course, experiments would detect a reduced number after 
accounting for efficiencies and acceptances. By itself, this 
number provides important insight into the driving force 
behind the disappearance of away-side jets. If the dis- 
appearance of away-side jets was largely an initial-state 
effect, e.g., 2 — ► 1 gluon fusion processes in the gluon 
saturation model [l2j, few extra particles would be cre- 
ated, especially at mid-rapidity. One would then expect 
to see a reduction in AN for Au+Au vs pp rather than 
the increase found here. 

If we assume a perfect detector, the momentum in- 
tegration described in Eqs. © and JSJ), resulted in 
AP X = —4.2 GeV/c, after triggering on leading hadrons 
in the p^ ls = 4—6 GeV/c range. This verifies momentum 
conservation in the procedure. As with the calculation 
of AN, most of this balancing momentum was found in 
the mid-rapidity cut. If the mechanism for jet quenching 
were an initial-state effect, one would expect less of the 
balancing momentum to be found with a rapidity near 
that of the triggering hadron. 

It would be straight-forward to expand this analysis to 
include a binning in rapidity as well as p x . One could 
then analyze the degree to which a jet spreads its energy 
and momentum in rapidity. Another possibility would be 
to bin with px or Et which would allow the determina- 
tion of the extra energy related to the trigger. The cal- 
culations shown here were performed with approximately 
4x 10 4 trigger particles. Two-dimensional binnings would 
require significantly higher statistics. We would certainly 
expect large-acceptance experiments such as STAR or 
PHENIX to be able to perform a comparable analysis 
as to what was demonstrated here if they were able to 
measure on the order of 10 5 events with trigger particles. 



Given sufficient statistics, the principal difficulty of 
these sort of analyses lies in defining a class of similar 
events for the purpose of background subtraction. For 
instance, if one uses multiplicity at mid-rapidity to de- 
fine the events, the quantity AN would always be zero. 
The means for determining centrality must be indepen- 
dent of the particles used in the analysis of p and po- 
Furthermore, it must be sufficiently accurate so that the 
observation of the triggering jet does not bias the event 
towards being more central. It might be quite possible to 
quantitatively estimate and correct for such a bias. For- 
tunately, these effects do not determine the extraction of 
AP X since the base distribution, po does not contribute 

We also investigated the degree to which our predic- 
tions are sensitive to the parameters used in the model. 
By reducing the medium-induced energy loss by a fac- 
tor of two, the net number of extra radiated particles 
was only reduced from 5.5 to about half that number. 
The energy loss is found to be rather insensitive to the 
parton elastic cross section. Once the medium becomes 
opaque, the behavior should be fairly insensitive to de- 
tails as thermalization should mask any details of the 
microscopic physics, and the linear dependence of the in- 
duced multiplicity with the energy loss parameters should 
dampen. Consequently, with a much higher initial gluon 
density, e.g., dN„/dv » 1000 is suggested in some gluon 
saturation models |25j , it is not clear to what degree the 
distribution of produced particles would be increased. 

In conclusion, we have shown that an analysis of low 
Pt observables can unveil the fate of lost jets at RHIC. 
A similar analysis of experimental data could provide 
irrefutable evidence for establishing final-state jet en- 
ergy loss as the explanation for jet suppression. A high- 
statistics multi-dimensional analysis would provide quan- 
titative insight into the dynamics and dissipation of jets 
in high-density matter. 
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